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We illustrate the magnetoelectric effect conditioned by the Katsura-Nagaosa-Balatsky (KNB) 
mechanism within the frames of exactly solvable spin-1/2 XY chains. Due to three-spin inter¬ 
actions which are present in our consideration, the magnetization (polarization) is influenced by 
the electric (magnetic) field even in the absence of the magnetic (electric) field. We also discuss 
a magnetoelectrocaloric effect examining the entropy changes under the isothermal varying of the 
magnetic or/and electric field. 
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I. INTRODUCTION 

Magnetically induced ferroelectricity is an important 
and highly interesting phenomena in condensed matter 
physics which attracts essential amount of attention both 
form experimentalists and theoreticianSfii^ In general, 
the phenomenon is a part of the so-called magnetoelec¬ 
tric effect, the intercoupling of magnetization and polar¬ 
ization in matter, which in its most prominent form can 
be briefly described as the magnetization dependence on 
the electric field and the polarization dependence on the 
magnetic field. The importance of the topic is not only 
limited to the highly nontrivial and rich physics, but it 
has very promising practical applications in future elec¬ 
tronic devices 41 ^ 

Among the various physical realizations of the magne- 
toelectric effect, the ferroelectricity of spin origin is par¬ 
ticularly important for spin-related electronics. There 
are several mechanisms of the coupling between local 
spins (magnetic moments) of the magnetic materials and 
the local polarization of the lattice celli^— The one, which 
we are going to exploit in the present study is called the 
Katsura-Nagaosa-Balatsky (KNB) mechanism,—!^ which 
is based on a so-called inverse Dzyaloshinskii-Moriya 
model or spin current model. This mechanism gives for 
the polarization (dipole moment) of the bond Pij con¬ 
necting two spins Si and the following expression; 

Py oc [eij X [sj X Sj]], (1.1) 

where e^- is the unit vector pointing form i-th site to j- 
th site. In terms of the spin current flowing between the 
magnetic sites 

jij oc [si X Sj], (1.2) 

the bond polarization in Eq. (HI) can be rewritten as 

Pij oc [eij X j,j]. (1.3) 


According to Eqs. (11.11) - (II. 3|) . a spiral spin order may 
generate a macroscopic polarization of electronic (spin) 
origin. 

At the present time a number of compounds is 
known which give evidence of magnetically-driven fer¬ 
roelectricity and magnetoelectric effect due to KNB 
mechanism.— The simplest quantum spin model ex¬ 
hibiting the spiral magnetic order is a spin-1/2 Ji- 
J 2 chain with ferromagnetic nearest-neighbor coupling 
(Ji) and antiferromagnetic next-nearest-neighbor cou¬ 
pling (J 2 ) which is considered to be more or less realistic 
model for such materials as LiCu 202 )^“— LiCuV 04 ^“— 
and CuCl 2 ^ to mention just few of them. 

In the present paper we are going to examine an exactly 
solvable model exhibiting a nontrivial magnetoelectric 
effect. The introduced quantum spin-chain model con¬ 
tains, in addition to the common two-spin interactions, 
the three-spin interactions and the bond polarization is 
of the spin origin according to KNB scenario (HD - dESD. 
The spin system under consideration is a variant of the 
famous spin-1/2 XY chaini^ (Suzuki model o^^d'^ ) which 
is solvable by the Jordan-Wigner fermionization. Such 
free-fermion spin models are quite popular as they pro¬ 
vide suitable playgrounds for the exact description of var¬ 
ious phenomena in strongly correlated systems— It is 
also in order to mention here similar recent studies on the 
exact treatment of the magnetism-driven ferroelectricity 
in quantum spin chains which, however, deal with the 
two-spin interactions only in the XXZ Heisenber g^^i^^ 
or compass modelsi^ The polarization for those mod¬ 
els is influenced by the magnetic field only in the pres¬ 
ence of the electric field, but the polarization is zero at 
zero electric field independently on the magnetic field. 
In other words, those models exhibit only the so-called 
trivial magnetoelectric effect. Contrarily to those pre¬ 
viously studied models, the ones considered below show 
the nontrivial magnetoelectric effect: The polarization is 
affected by the magnetic field even at zero electric field. 

The rest of the paper is organized as follows. In 
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Sec.ini we discuss the spin-chain models of magnetism- 
driven ferroelectricity which can be treated exactly by 
the Jordan-Wigner transformation to free fermions. Ex¬ 
act solutions are presented in Sec. uni The essential fea¬ 
ture of the models at hand is the three-spin interactions. 
For the sake of simplicity we distinguish two types of 
the three-spin interactions, the XZY — YZX ones and 
the XZX + YZY ones, which are considered separately 
in Sec. 113 and Sec. El respectively. Our main focus 
is on the ground-state magnetization and polarization 
(zero-temperature properties) as well as on the magneto¬ 
electrocaloric effect (finite-temperature properties). We 
summarize our findings and sketch perspectives for fur¬ 
ther work in Sec. ED Details of some calculations are 
presented in two appendices. 


II. SPIN-1/2 XY CHAINS AUGMENTED BY 
KNB MECHANISM 

In this paper, we consider N spins 1/2 placed on a 
one-dimensional linear-chain lattice. We start with the 
“bare” Hamiltonian of the spin system 


Ho = H,j + He+Hk + H^, (2.1) 

which contains, in general, the ordinary two-spin 
isotropic XY interactions, 

Hj = jY.i«+i + ^n4+i) (2-2) 


(we may set J > 0 without sacrificing generality), the 
three-spin interactions of two typeSj ^^d^ 

He = eY, (4<+i4+2 - 4 <+ 1 <+ 2 ) , 


Hk = kJ2 («+i <+2 + , (2-3) 


as well as the Zeeman interaction of the spins with a 
z-aligned external magnetic field 'H, 

Hz = -nJ2<- (2-4) 


Boundary conditions are not considered explicitly for 
thermodynamically large systems iV —> oo. 

In our study, we assume that the spin model arises from 
a more fundamental electronic model for which the KNB 
scenario^ holds. This means that the electric polarization 
of the bond between the neighboring sites n and n -|- 1 to 
be denoted by p^.n-i-i is determined by the spin current 
in,n+i according to Eqs. (II.ip - (II.3p . In our model, the 
chain runs along the x-direction in the real space, i.e., 
Gij = (1,0,0) in Eqs. (11.11) . (II.3p . Furthermore, the x- 
, y-, and z-axis in the real space and in the spin space 
coincide. Then p„,„+i oc (0,„+i), i.e., the 
bond polarization has zero x-component. 


The spin-current operator satisfies the lattice version 
of the continuity equation, ds“/dt = —divj“, a = x, y, z, 
where the discrete divergence operator acts in the real 
space. Consider at first the isotropic XY interactions 
only [i.e., E = A = 0 in Eqs. (12.11) - (12.41) ]. For the 
z-component of the current through the bond we find: 

^ = -i [<, Hj +i/z] = - (j:+. - .) = -divj:, 

(<4+1-44+0 ■(2-5) 

If the spin-1/2 isotropic XY chain comes from an un¬ 
derlying electronic model for which KNB mechanism 
works, then each bond has the polarization Pn,n-i-i oc 
(0> -jn,n+i,jl,n+i) ^ith < „+i = <+1 given in Eq. ([21]) 
which may manifest itself in the presence of electric field. 
Let £ be a y-aligned external uniform electric field [i.e., 
the electric field vector is (0,£’, 0)]. Then the Hamilto¬ 
nian of the model Hq has to be supplemented with the 
term -ST.nPl.,n+i ^£Y.njn,n+i and we arrive at^* 

h = jY. ^ - 44+l)] 

n 

-71 ^4. (2.6) 


One may calculate the magnetic moment XP = 
and the polarization oc —J J2n(^n^n+i ~ ^n^n+i) 
for spin model (12.61) using the standard Jordan-Wigner 
fermionization method,— see Ref. [2^ and Sec. IbD 
Next we consider the model with the three-spin in¬ 
teractions, see Eq. (|2.3I) . Placing the system in a uni¬ 
form external electric field (0,£’,0) we have to add to 
the Hamiltonian Hq the term —SP^. Following previ¬ 
ous studies on the magnetoelectric effect in the spin-1/2 
J1-J2 chainj^^— we take for the required polarization 
operator the following expression: P^ = J2nPn n+i 
- 'Enjn,n +1 ^ith < „+i = <+1 given in Eq. (EH). Then 
the Hamiltonian of the spin system becomes 

h = jJ2 [44+1 + 44+1 + f (44+1 - 44+i)] 

n 

+e^(44+is](+2-44+i4+2) 

n 

(sX+i4+2+ 44+i4+2) 

n 

-J^^4.(2.7) 


As previously, the magnetic moment = 3n('®n) 
and the polarization P^ oc Z[n(®n'Sn+i “ ^n^n+i) for 
spin model (12.71) can be calculated rigorously within the 
frames of the Jordan-Wigner fermionization approach^ 
see Sec. ED It should be also noted here, that models sim¬ 
ilar to the one given in Eq. (1^ have been considered 
recently [compare Eq. (12.71) to Eq. (2.1) of Ref. or to 
Eq. (2.2) of Ref. however, they contain the two-spin 



3 


interactions only. As it is shown below, the three-spin 
interactions lead to new features of the magnetoelectric 
effect. 

Interestingly, for the spin model defined in Eqs. (EH) 
- (12.41) we may also consider the three-spin-interaction 
contribution to the spin current in KNB formula for the 
bond polarization (11.31) (see Ref. and remain to face 
an exactly solvable spin model. Indeed, using the conti¬ 
nuity equation, for the required z-component of the spin- 
current operator we find: 

^ = -i [<, ffj, + Hjc] = = -divj^, 

j^+l = -2E (s“4+iS“+2 + 44-Hi4-H2) 
+2K - 4<+i<+2) ■ (2.8) 

Therefore in the presence of the electric field (0,£1, 0) we 
face the Hamiltonian similar to the one given in Eq. (123) 
but with 


E^E + 2K£, K^K-2E£. 


(2.9) 


Clearly, all averages are calculated now with this new 
Hamiltonian (12.7L (12.91) . The magnetic moment is given 
by whereas the polarization now is 

the sum of the two terms, oc —JJ2ni^n^n+i ~ 
s((s^+l) and py oc J 2 n(‘^^(SnK+lSn +2 + 'SJK-Hl 4-^2) “ 
2 A:(s“s^+i4+2 - sls^+iS^+ 2 )), see Eq. ([13. Again 
these quantities can be calculated rigorously using the 
Jordan-Wigner transformation to fermionSfi^ see Sec. lHIl 

In summary, basing on the KNB scenario for 
magnetism-driven ferroelectricity (fTTl) - (fOl) and choos¬ 
ing specific exchange interactions, lattice geometry, and 
direction of external fields we arrive at simple spin 
models which are expected to mimic basic features of 
some multiferroics. In contrast to the considered pre¬ 
viously similar exactly-solvable models of magnetism- 
driven ferroelectricity,—^— we have included into con¬ 
sideration also the three-spin interactions. More real¬ 
istic description of multiferroics of spin origin leads to 
more sophisticated spin models (see, e.g.. Refs. IsTI-l^ 
and their further analysis involves approximations. The 
merit of the introduced models is their exact solvability. 
In the following sections we examine rigorously some as¬ 
pects of the magnetoelectric effect in spin models defined 
in Eq. (12.61) . Eq. (12.7L and Eqs. (12.71) . (12.91) . 

To close the section, we make a remark concerning the 
three-spin interactions which are present in Eq. ism . 
Obviously, the spin model with two-spin interactions only 
(1^ describes the influence of £ on or of T-L on 
py. However, from symmetry arguments it is clear that 
= 0 if "H = 0 independently on f or = 0 if f = 0 
independently on "H, see the paragraph after Eq. (1^ in 
sec.nni In contrast, in the presence of the three-site in¬ 
teractions (12.31) these symmetry arguments do not work 
any more and one may expect a nontrivial magnetoelec¬ 
tric effect when IvP is influenced by £ even for P = 0 or 
py is influenced by % even for f = 0. Rigorous calcula¬ 
tions reported in Secs. IIVI and IVl support these expecta¬ 
tions. 


III. EXACT SOLUTIONS 


Spin-1/2 models given in Eqs. (12.6|) . (12.7|) . and (12.7L 
(1^ are exactly solvable via the Jordan-Wigner transfor¬ 
mation to fermionSfi^ see Appendix lAl In fermionic pic¬ 
ture we face noninteracting spinless fermions with known 
energies and thus many statistical mechanical calcu¬ 
lations can be easily carried out. For the Helmholtz free 
energy per site we find^^i^l 

f(T,n,£) = -^ £ d/tln(2cosh|^) , (3.1) 

where for model (123 

Ck = —H + J cos K + J£ sin k 
E K 

— — sin(2 K)—^cos(2k). (3.2) 

In the case of model (12.6|) one has to put A = AT = 0 in 
Eq. ([33, whereas in the case of model (123, dlH) one 
has to make the replacement ()2.9D in Eq. (13.2p . 

In what follows we are interested in the magnetization 
and the polarization for the spin models at hand. For 
the magnetization per site we have 



0 1- (3.3) 

1+ tT 

For the polarization per site for spin model (12.71) we have 
p = py (X - 44 - hi ) 

n 

J r 

= —— / dAtsinKu^. (3-4) 


While calculating the polarization per site for spin model 
(1271) . (123 we have to use the energy spectrum (IQ) . 
(1^ and in addition to the contribution py given in 
Eq. (13.41) to calculate the second term p^. 


p = oc 


1 




^ ^ (2A(s^s^+iS^+2 + 4<-hi4-h2) 

n 

-2K(«+i4+2 - 4<+i4+2)) 

1 r 

/ dK [Ecos(2k) — Arsin(2K)] tIk. (3.5) 
27r 


Furthermore, we are interested in the entropy (per site) 
for the spin models at hand. It immediately follows from 
Eq. (13.1|) through the relation s = —df/dT, 


s = — / dK 
27r 


_ln(2coshg;)-|^tanh^]. (3.6) 


As emphasized in the end of the previous section, the 
introduced three-spin interactions are essential for ap¬ 
pearance of the nontrivial magnetoelectric effect. This is 
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clearly seen from formulas (lO) - dH in the fermionic 
picture. Obviously, for T = 0 we have = 0 if e„ > 0 
and = 1 if e„ < 0. Therefore, for E = K = 0 we 
have d/CTT-K = tt if 'H = 0 independently on £, and 
hence the ground-state magnetization (j3.3|) is zero. Fur¬ 
thermore, for E = 0 we have if £ = 0 inde¬ 

pendently on 'H, and hence the ground-state polarization 
(EH) is zero. In the presence of the three-site interactions 
the used symmetries of (13.21) may be broken and the 
obtained conclusions do not hold any more. 

For the sake of simplicity in the following analysis 
we distinguish two different types of three-spin interac¬ 
tions, XZY — YZX type, i.e., K = Q (Sec. IIVI) and 
XZX + YZY type, i.e., E = 0 (Sec. 0. From previous 
studies^^— we know that for f = 0 both models 
[see Eqs. (12.11) - (12.41) ] exhibit three phases in the ground 
state: High-held ferromagnetic phase and, depending on 
the relation between two- and three-spin interactions, two 
different spin-liquid phases. In the fermionic picture dif¬ 
ferent phases correspond to different Fermi-surface topol¬ 
ogy (different number of Fermi points) for fermions. Fur¬ 
thermore, these models are known to exhibit a nonzero 
magnetization in zero magnetic heldJ^d^ Clearly, for 
£ ^ 0 both types of the three-spin interactions are 
present in the analysis of the magnetoelectric effect. 
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IV. THREE-SPIN INTERACTIONS OF 
XZY - YZX TYPE 

In this section we consider the case of the XZY — 
YZX three-spin interaction, i.e., if 0, AT = 0 (J = 1). 
From the ground-state phase diagram of the model (12.71) 
with £ = Q (see, e.g.. Ref. [2^ we know that it makes 
sense to distinguish three representative values of if, for 
example, E = 0.5, if = 1, and E = 2. For these cases 
we calculate the ground-state phase diagram in the plane 
'H-£ which indicates phases having different numbers of 
Fermi points [see Fig. [Tjfor model (12.71) with E = 2], and 
the ground-state^ magnetization and polarization [see 
Fig. [2| for model (12.71) with E = 0.5 and E = 2: m{'H) 
(bold solid), p(ii) (bold dashed) at £ = 0 and m{£) (bold 
solid), p{£) (bold dashed) at ii = 0]. Furthermore, we 
calculate the low-temperature entropy in the plane 'H~£ 
[see Fig. Ofor model (12.71) with if = 2 at T = 0.09]. 

Let us discuss the obtained results. The physics of 
the model along the line £1 = 0 in the plane 'H-£ is well 
understood While TL increases from —oo to oo the 

system passes from the ferromagnetic phase to the ferro¬ 
magnetic phase over the spin-liquid I phase if |if| < 1 or 
over the spin-liquid I phase, the spin-liquid II phase, and 
again the spin-liquid I phase if |if| > 1. The sequences of 
phases for the latter case with E = 2 can be followed in 
Fig. m along the line £ = {). In the presence of the electric 
field £ ^ Q, the system exhibits the same phases sepa¬ 
rated by quantum phase transition lines and two quan¬ 
tum triple point. Ground-state dependences m{'H) (solid 
lines) and p{TL) (dashed lines) at f = 0 for A = 0.5 and 


FIG. 1: (Color online) Magnetization (upper panel) and po¬ 
larization (lower panel) at very low temperature (T = 0.005) 
along with the ground-state phase diagram (black solid lines 
separate different phases) of model (12.711 with J = 1, E = 2, 
K = 0. 




FIG. 2: (Golor online) Dependencies of magnetization (bold 
solid) and polarization (bold dashed) at T = 0.005 on mag¬ 
netic field at £1 = 0 (first row) and on electric field at 'H = 0 
(second row) for model (12.711 with J — 1, E = 0.5 (left col¬ 
umn) or E — 2 (right column ), K = 0. By thin lines we show 
magnetization (thin solid) and polarization (thin dashed) for 
model (12.71) . (12.91) with J = 1, E = 0.5, 2, K = 0. 
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FIG. 3: (Color online) Towards magnetoelectrocaloric effect 
for model (12.711 with J = 1, i? = 2, = 0: Entropy de¬ 

pendence on magnetic and electric fields at low temperature 
T = 0.09 along with the ground-state phase diagram (black 
solid lines separate different phases). 


E = 2 are shown in the first row in Fig. [2j These depen¬ 
dences are trivial only in the ferromagnetic phase, when 
m{'H) = ±1/2 and p('H) = 0. In the both spin-liquid 
phases, not only the magnetization is influenced by the 
magnetic field but also the polarization even for £ = 0 
is effected by the magnetic field. Furthermore, from the 
ground-state dependences m[£) and p{£) at "H = 0 (the 
second row in Fig. [2|) one can see that in the both spin- 
liquid phases not only p depends on £ but also m even 
for 'H = 0 is governed by £. 

For the model at hand we can easily obtain rigorously 
further details about the behavior of relevant quanti¬ 
ties. Consider, e.g., the ground-state polarization p('H) 
slightly below the saturation field 'Hsa.t > 0 for model 
(12.71) with J = 1, E = 0.5 at £1 = 0, see the corre¬ 
sponding panel in Fig. [2] Clearly, the saturation field 
"Hsat is defined by the fermion energy spectrum (lO) . 
namely, —'Hsa.t ± cosk* — 0.25 sin(2A{*) = 0, where k* is 
determined from the equation deK/dK\K,=K,* = — sinK* — 
0.5 cos(2k*) = 0, i.e., k* Ri —0.374 734. Therefore 
'Hsat ~ 1.100 917. Assume further 'H = 'Hsat — 5'H, where 
5'H > 0 is a small quantity. According to Eq. (13.4L 
p = (coskJ— cosK2)/(27r), where k\ r; k* —1.113 OlhV^'H 
and /t 2 ~ K* ± 1.113 915-\/^. As a result, at 77 = 
'Hsat — S'H we have p('H) r; —0.129 782V77sat — 77, that 
is, the ground-state polarization emerges with the criti¬ 
cal exponent 1/2 as the control parameter 77 passes the 
critical value 77sat- 

The ground-state dependences of m and p ± p on 77 
and £ for model (12.7L (12.911 are shown in Fig. [2] by thin 
lines. Within the fermionic picture (13.21) . (13.41) . (13.51) it 
can be proved (Appendix I bI) that for this model p±p = 0 
if f = 0 independently on 77, see thin dashed lines in the 
two panels from the first row in Fig. [2l For nonzero 
however, p±p is influenced by 77. Clearly, to at f = 0 for 
the model (12.71) and the model (12.71) . (12.91) is the same. 

The ground-state phase diagram also manifests itself 


in the dependence of the entropy s (13.61) on 77 and £ 
at low temperatures. The low-temperature entropy ex¬ 
hibits well pronounced maxima along the quantum phase 
transition lines and around the quantum triple points. 
These maxima become more sharper as the tempera¬ 
ture decreases. If the system is placed in a thermo¬ 
stat with the temperature T, AQ = T AS” with AS” = 
S'(772,£’2) — <S'(77i, fi) is the heat the system takes in (if 
AS” > 0) or gives out (if IS.S < 0) under the change of the 
fields from the values 77i, f i to the values 772 , £^2 • Clearly, 
the system at hand exhibits a magnetoelectrocaloric ef¬ 
fect, i.e., can be used for cooling/heating under a change 
of external fields. The magnetoelectrocaloric effect is 
most pronounced at low temperatures around the quan¬ 
tum phase transition lines and around the quantum triple 
points, see Fig. [31 


V. THREE-SPIN INTERACTIONS OF 
XZX + YZY TYPE 

We pass to the case of the XZX ± YZY three-spin 
interaction, i.e., E = 0, K ^ 0 {J = 1). From the 
ground-state phase diagram of the model (123) with £ = 
0 (see, e.g., Ref. [13) we know that it makes sense to 
distinguish three representative values of K, for example, 
K = 0.25, K = 0.5, and K = 1.7. For these cases we 
calculate the ground-state phase diagram in the plane 
77which indicates phases having different numbers of 
Fermi points [see Fig. [Ufor model (12.71) with K = 1.7], 
and the ground-state^ magnetization and polarization 
[see Fig. [S]for model (12.71) with K = 0.5 and K = 1.7: 
to( 77) (bold solid), p(77) (bold dashed) at £i = 0 and m{£) 
(bold solid), p{£) (bold dashed) at 77 = 0]. Furthermore, 
we calculate the low-temperature entropy in the plane 77~ 
£ [see Fig.[6]for model (12.71) with 7F = 1.7 at T = 0.09]. 

Again the properties of the model along the line £ = 0 
are well known<22i21 As 77 varies from —00 to 00 , the 
system is driven by varying 77 from the ferromagnetic 
phase to the ferromagnetic phase through the spin-liquid 
I phase if \K\ < 1/2 or through the spin-liquid I and 
spin-liquid II phases (through the spin-liquid II and spin- 
liquid I phases) if K > 1/2 (if K < —1/2). Ground- 
state dependences of to and p (and p ± p) on fields show 
nontrivial features only outside the ferromagnetic phase. 
However, p(77) = 0 along the line £ = 0 that is obviously 
traced back to the symmetry = e_K, see the first row 
in Fig. [5l Breaking this symmetry by switching on £ 
immediately results in influence of 77 on p, see the second 
row in Fig. [5| On the other hand, to is influenced by £ 
even at 77 = 0, see the third row in Fig. [5] Again the 
considered model exhibits a magnetoelectrocaloric effect, 
which is most pronounced at low temperatures around 
characteristic lines of the ground-state phase diagram, 
see Fig. [5] 
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FIG. 4: (Color online) Magnetization (upper panel) and po¬ 
larization (lower panel) at very low temperature {T = 0.005) 
along with the ground-state phase diagram (black solid lines 
separate different phases) of model (12.71) with J = 1, = 0, 

K = 1.7. 
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FIG. 5: (Color online) Dependencies of magnetization (bold 
solid) and polarization (bold dashed) a.t T — 0.005 on mag¬ 
netic field ai £ = 0 (first row) and £ = 0.7 (second row) and 
on electric field at 'H = 0 (third row) for model (12.71) with 
J = 1, E = 0, K = 0.5 (left column) or K — 1.7 (right col¬ 
umn). By thin lines we show magnetization (thin solid) and 
polarization (thin dashed) for model (12.71) . (12.91) with J = 1, 
E^O, K = 0.5, 1.7. 



VI. CONCLUSIONS 


We have considered simple but nontrivial models of a 
multiferroic of spin origin. The main worth of the mod¬ 
els is their exact solvability: All relevant quantities can 
be calculated rigorously and examined in detail. These 
studies may serve as a benchmark for more realistic cases 
which are not exactly solvable. 

In contrast to free-fermion models studied earlier 
we include in the model the three-spin interactions of 
XZY — YZX and XZX + YZY types. Due to these 
interactions the magnetoelectric effect becomes espe¬ 
cially interesting: Magnetization (polarization) can be 
induced and governed solely by electric (magnetic) field. 
The considered models show magnetoelectrocaloric ef¬ 
fect, i.e., isothermally (adiabatically) varying fields no¬ 
ticeably change the entropy (temperature). The effect 
is most pronounced at low temperatures around pecu¬ 
liarities (quantum phase transition lines and quantum 
triple points) on the ground-state phase diagram which 
is rather rich in the presence of the three-spin interac¬ 
tions. 

It should be stressed that some characteristic features 



‘X 


FIG. 6: (Golor online) Towards magnetoelectrocaloric effect 
for model (12.71) with J = 1, E = Q, K = 1.7: Entropy de¬ 
pendence on magnetic and electric fields at low temperature 
T = 0.09 along with the ground-state phase diagram (black 
solid lines separate different phases). 
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of the reported dependences for the polarization or the 
magnetization on fields (e.g., emergence after passing a 
threshold value, cusps, abrupt changes, etc.) can be seen 
in experimentally measured data, see Ref. and refer¬ 
ences therein. On the other hand, the minimal model 
to describe such spin-chain multiferroics as LiCu 202 or 
LiCuV 04 is the spin-1/2 anisotropic Heisenberg model 
with the Hamiltonian 

= («+1 + + A < 4 + i ) 

n 

+ J2 (SnSn+2 + SnSn+2 + As* S^_|_2) 

- 774 ], (6.1) 


to get for the Hamiltonian (12.7p a bilinear Fermi-form: 


H 


E 


J 0. lEK I 

2 ^n^n+l ^ C^Cn+2 + H.C. 


-U 1 cicn - - 


.(A2) 

This Hamiltonian can be brought into the diagonal form 
after the Fourier transformation 


E E (A 3 ) 


with ferromagnetic Ji < 0, antiferromagnetic J 2 > 0, 
and small easy-plane anisotropy A < 1, see, e.g.. Ref. H 
One way to examine this model is to apply the Jordan- 
Wigner fermionization method. However, in contrast to 
the model given in Eqs. (ED) - (EH), one arrives for 
model EU at interacting fermions and further treatment 
becomes approximate. 

Finally, the considered models hold promise as a core 
system permitting to examine some other aspects of mul¬ 
tiferroics, e.g., the dynamical magnetoelectric effecti^ 
The work in this direction is in progress. 
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Appendix A: Jordan-Wigner transformation and 
thermodynamic functions 

Let us consider statistical mechanical calculations for 
the spin-I/2 chain model given in Eq. (12.71) . First we 
introduce the operators = sj ± is^. Then we use the 
Jordan-Wigner transformation to spinless fermions, 

c\=st, c] = (-2sD..-(-2s|_i)s+,j=2,...,iV, 

Cl = S", Cj = (-2sJ)... (-2s|_i)sJ, j = 2,.. .,iV/Al) 


K = 2TTm/N, m = — A/2,..., A/2 — 1 (we assume with¬ 
out loss of generality that A is even), which yields 

K ^ ^ 

with Ck given in Eq. (I3.2|l . 

The partition function of A Fermi oscillators (1X41) can 
be easily calculated, Z{T,'H,£,N) = ]/[^ 2cosh[eK/(2T)]. 
It yields the Helmholtz free energy per site given in 
Eq. (13.11) . To get the magnetization (13.31) and the po¬ 
larization (13.41) one may simply take the corresponding 
derivatives, i.e., m = —df/d'H and p = —dfldS. 


Appendix B: Ground-state polarization for model 

(Em . (Em at f = 0 

Consider the ground-state polarization p+p for model 
(12.71) . (12.91) . According to Eqs. (13.41) . (13.51) we have to 
calculate 


— — / ds; [ J sin K J-A cos(2k) — A sin(2K)]. (Bl) 

27r Jp 

Here F is the domain within the interval [—tt, tt] where 
Ck = —H+J cos K—(A/2) sin(2K) —(A/2) cos(2k) < 0 (we 
consider the case £ = 0). The antiderivate in Eq. (IBII) 
is — J cosK -I- (A/2) sin(2K) -|- (Kj^') cos(2«:) = —— A 
and therefore the resulting integral is zero since 6^ = 0 
at the boundaries of the domain F. Moreover, p + p — 0 
for model (12.71) . (j2.9l) at £ = 0 for nonzero temperatures 
too. 


^ M. Fiebig, Journal of Physics D: Applied Physics 38, R123 ^ Y. Tokura, S. Seki, and N. Nagaosa, Rep. Prog. Phys. 77, 

(2005). 076501 (2014). 







^ S.-W. Cheong and M. Mostovoy, Nature Materials 6, 13 
(2007). 

^ Y. Tokura and S. Seki, Advanced Materials 22, 1554 

( 2010 ). 

® H. Katsura, N. Nagaosa, and A. V. Balatsky, Pliys. Rev. 
Lett. 95, 057205 (2005). 

® C. Jia, S. Onoda, N. Nagaosa, and J. H. Han, Phys. Rev. 
B 74, 224444 (2006). 

^ S. Park, Y. J. Choi, C. L. Zang, and S.-W. Cheong, Phys. 
Rev. Lett. 98, 057601 (2007). 

® S. Seki, Y. Yamasaki, M. Soda, M. Matsuura, K. Hirota, 
and Y. Tokura, Phys. Rev. Lett. 100, 127201 (2008). 

® A. A. Bush, V. N. Glazkov, M. Hagiwara, T. Kashiwagi, 
S. Kimura, K. Omura, L. A. Prozorova, L. E. Svistov, 
A. M. Vasiliev, and A. Zheludev, Phys. Rev. B 85, 054421 
( 2012 ). 

“ Y. Qi and A. Du, Phys. Lett. A 378, 1417 (2014). 

Y. Naito, K. Sato, Y. Yasui, Y. Kobayashi, and M. Sato, 

J. Phys. Soc. Jpn. 76, 023708 (2007). 

Y. Yasui, Y. Naito, K. Sato, T. Moyoshi, M. Sato, and 

K. Kakurai, J. Phys. Soc. Jpn. 77, 023712 (2008). 

F. Schrettle, S. Krohns, P. Lunkenheimer, J. Hemberger, 

N. Biittgen, H.-A. Krug von Nidda, A. V. Prokofiev, and 
A. Loidl, Phys. Rev. B 77, 144101 (2008). 

S. Seki, T. Kurumaji, S. Ishiwata, H. Matsui, H. Mu- 
rakawa, Y. Tokunaga, Y. Kaneko, T. Hasegawa, and 
Y. Tokura, Phys. Rev. B 82, 064424 (2010). 

E. Lieb, T. Schultz, and D. Mattis, Annals of Physics 16, 
407 (1961). 

M. Suzuki, Phys. Lett. A 34, 94 (1971). 

M. Suzuki, Prog. Theor. Phys. 46, 1337 (1971). 

D. Gottlieb and J. Rossler, Phys. Rev. B 60, 9232 (1999). 

O. Derzhko, J. Richter, and V. Derzhko, Annalen der 
Physik 8, SI-49 (1999) arXiv:cond-mat/9908425 . 

1. Titvinidze and G. 1. Japaridze, Eur. Phys. J. B 32, 383 


(2003). 

P. Lou, W.-C. Wu, and M.-G. Chang, Phys. Rev. B 70, 
064405 (2004). 

A. A. Zvyagin, Phys. Rev. B 72, 064419 (2005). 

T. Krokhmalskii, O. Derzhko, J. Stolze, and 

T. Verkholyak, Phys. Rev. B 77, 174404 (2008). 

O. Derzhko, T. Krokhmalskii, J. Stolze, and 

T. Verkholyak, Phys. Rev. B 79, 094410 (2009). 

V. Derzhko, O. Derzhko, and J. Richter, Phys. Rev. B 83, 
174428 (2011). 

F. G. Ribeiro, J. P. de Lima, and L. L. Gongalves, J. Magn. 
Magn. Mater. 323, 39 (2011). 

M. Topilko, T. Krokhmalskii, O. Derzhko, and 

V. Ohanyan, Eur. Phys. J. B 85, 27 (2012). 

M. Brockmann, A. Kliimper, and V. Ohanyan, Phys. Rev. 
B 87, 054407 (2013). 

P. Thakur and P. Durganandini, AIP Gonference Proceed¬ 
ings 1665, 130051 (2015). 

W. -L. You, G.-H. Liu, P. Horsch, and A. M. Oles, Phys. 
Rev. B 90, 094413 (2014). 

C. Jia and J. Berakdar, Europhys. Lett. 85, 57004 (2009). 
C. Jia and J. Berakdar, Journal of Superconductivity and 
Novel Magnetism 25, 2679 (2012) arXiv:1101.2067 . 

M. Azimi, L. Chotorlishvili, S. K. Mishra, T. Vekua, 
W. Hiibner, and J. Berakdar, New J. Phys. 16, 060318 
(2014). 

J. Sirker, Phys. Rev. B 81, 014419 (2010). 

Although the results reported in Figs. [2] and [5] refer to a 
very small but nonzero temperature (T = 0.005) they can 
not be distinguished from the results at T = 0 in the scale 
of the figures. 

36 Y, Tokura and N. Kida, Phil. Trans. R. Soc. A 369, 3679 

( 2011 ). 



